Stacking-order dependent transport properties of trilayer graphene 
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We report markedly different transport properties of ABA- and ABC-stacked trilayer graphenes. 
Our experiments in double-gated trilayer devices provide evidence that a perpendicular electric field 
opens an energy gap in the ABC trilayer, while it causes the increase of a band overlap in the ABA 
trilayer. In a perpendicular magnetic field, the ABA trilayer develops quantum Hall plateaus at 
filling factors of v — 2, 4, 6... with a step of Ais = 2, whereas the inversion symmetric ABC trilayer 
exhibits plateaus at v — 6 and 10 with 4- fold spin and valley degeneracy. 



PACS numbers: 73.43.-f, 72.80.Vp, 73.63.-b 

The unique chiral nature of low-energy quasiparticles 
in graphene, characterized by a Berry phase J it with 
linear and parabolic dispersion for monolayer (J = 1) 
and bilayer ( J = 2) graphene respectively, results in un- 
usual quantum Hall effects (QHE)PH^The Landau-level 
(LL) energy in a perpendicular magnetic field £?, given 
by E n oc y/Bn for monolayer and E n oc B^/n(n — 1) for 
bilayer graphene, shows J-fold degenerate LLs at zero 
energy, with integer n being the LL index. The existence 
of J-fold degenerate zero-energy LLs, combined with 4- 
fold spin and valley degeneracy in each LL, explains the 
unusual sequence of quantum Hall states observed at fill- 
ing factor sequences v — ±2, ±6, ±10... for monolayer}^ 
and v — ±4, ±8, ±12... for bilayer graphenePThe bilayer 
graphene is further distinguished from the gapless mono- 
layer by a tunable energy gap, induced by breaking the 
inversion symmetry of the two layers in a perpendicular 
electric neld.EHSl 

In few-layer graphene (FLG), the stacking order of- 
fers an extra degree of freedom. Indeed, the electronic 
structure and the Landau level spectrum differ signif- 
icantly depending on the stacking order in FLGP^EQ 
For instance, Bernal (ABA) -stacked trilayer exhibits an 
electric-field tunable band overlap while rhombo- 
hedral (ABC)-sta cked tr ilayer is predicted to present a 
tunable band gap! 7 * 10 * 1 -^ To date, no direct evidence of 
the electric-field and stacking-order dependent transport 
properties has yet been reported in double-gated devices. 
In the simplest tight-binding model that includes only the 
nearest intra- and inter-layer hopping parameters 70 and 
71 (Fig. 1(b)), the Landau level spectrum of the ABA tri- 
layer can be viewed as a superposition of V^-dependent 
monolayer- like LLs and Independent bilayer-like LLs 
(Fig. 1(c)). On the other hand, LLs of the ABC trilayer 
(Fig. 1(d)) are given by E n ex B z < 2 ^n{n - l)(n - 2) 
with Berry's phase 3tt. 7 8 Despite the substantial differ- 
ence in the LL spectrum, 3-fold degenerate zero-energy 
LLs with 4-fold spin and valley degeneracy are expected 
to result in QHE plateaus at filling factor sequences 
v = ±6, ±10, ±14... for the trilayer graphene indepen- 



dently of the stacking order E^ 14 * 15 * However, the lack of 
inversion symmetry in ABA trilayer may lead to a bro- 
ken valley degeneracy, while the valley degeneracy of LLs 
is guaranteed in the inversion-symmetric ABC trilayer.^ 
Here, we report stacking-dependent transport properties 
of double-gated trilayer graphene, combined with Raman 
spectroscopy. We show that the effects of applied elec- 
tric and magnetic fields on the ABC-stacked trilayers are 
strikingly different from those on the ABA-stacked tri- 
layers. 

Our experiments have been perfor med on single- and 
double-gated trilayer graphene devices^^ prepared by 
exfoliating graphite on Si/Si02 substrates. The heavily 
conductive Si was used as a back gate and the thickness 
of the Si02 layer was 285 nm. Ti/Au electrodes and top 
gates (Si02 /Ti/Au) were defined on top of the graphene 
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FIG. 1: (Color online) (a) Optical images of a typical device 
before and after top-gate deposition, (b) The structure of 
ABA-stacked trilayer graphene. (c) Landau levels of ABA- 
stacked trilayer graphene shown up to n — 10, calculated with 
70 = 3 eV and 71 = 400 meV. (d) Landau levels of ABC- 
stacked trilayer graphene shown up to n — 10. 
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flakes by electron-beam lithography (Fig. 1(a)). The 
dc magneto-transport properties were studied at liquid- 
Helium temperatures in pulsed perpendicular magnetic 
fields of up to 50 T. The magnetic-field pulse resulted 
from the discharge of a large capacitor bank with a ca- 
pacitance of 30 mF and a voltage up to 20 kV and lasted 
typically ~ 500 ms. 

All our investigated graphene samples have been reli- 
ably identified as trilayers and their stacking order has 
been determined by means of Raman spectroscopy. We 
used an excitation laser with a wavelength of 532 nm and 
a spot size of 1.5 fim in diameter. The Raman spectra of 
mechanically exfoliated graphene shows two peaks: the 
G band and the 2D (GM band at respectively 1580 cm -1 
and 2700 cm -1 (Fig. |2[a)). The G band is due to the 
first-order Raman scattering by the double degenerate 
E2 g phonon mode at the Brillouin zone center, while the 
2D band originates from a second-order process, involv- 
ing two intervalley optical phonons near the boundary of 
the Brillouin zoneP^ The peak at 520 cm -1 (labeled as 
Si) is due to the first-order Raman scattering by optical 
phonons of the Si substrate. 

A reliable approach to count the number of layers (TV) 
of FLG deposited on Si/Si02 substrates is based on the 
ratios of the intensities of the G peak and the Si peak, 
Ig/Isi- 19 As shown in Fig. [2^a), for a flake containing 
up to 7 layers, the intensities of the G and Si peaks 
clearly change with N. We find that Ig/Isi increases 
monotonically and discretely with N due to an increase 
of the intensity of the G peak and a decrease of the in- 
tensity of the Si peak. Our findings are in agreement 
with recent observations, which attribute this behavior 
to enhanced absorption and Raman scattering of light 
by thicker graphene layers. 19 In Figure (2^b), we show the 
Raman spectra of all trilayers investigated here. These 
samples have G and 2D peaks of similar intensities and 
their Iq /Isi is consistent with the typical values found for 
trilayer graphene. 

An accurate determination of N for up to three lay- 
ers is also possible from the 2D peak since its shape 
and position evolves with N (see Fig. [2ja)). The 2D 
band is affected by the band structure of the material 
since it arises from a double-resonance process involv- 
ing transitions among various electronic states. As tri- 
layer graphene has three valence and three conduction 
bands, up to 15 electronic transitions can contribute to 
the 2D band. 20 However, many of these different pro- 
cesses have very close energy separations and experimen- 
tally it is found that the minimum number of Lorentzian 
functions necessary t o cor rectly fit the 2D mode of tri- 
layer graphene is six. 20 21 Consistently, Fig. 2(c) shows 
that for all the investigated trilayer graphene samples a 
good fit can be achieved with 6 Lorentzian functions. 

Having determined the number of layers, we now ad- 
dress the stacking order in our trilayer samples. It has 
been recently demonstrated that an accurate and efficient 
method to characterize stacking order in FLG is based 
on the distinctive features of the Raman 2D peakP^ We 
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FIG. 2: (a) Raman spectra for graphene samples with dif- 
ferent number of layers. The inset shows the optical micro- 
scope picture of the flake containing up to 7 layers used for 
these measurements. The labels 1L to 7L indicate the num- 
ber of layers, (b) Raman spectra of trilayer samples discussed 
in this work, (c) The 2D Raman band of graphene trilayers 
with ABC stacking (TGI and TG2) and ABA stacking (TG3). 
The red lines are fits by 6 Lorentzian functions and the lines 
of other colors are the Lorentzian components of the fits. 



find that TGI and TG2 show a more asymmetric 2D 
peak than TG3, consistent with the reported differences 
between ABC and ABA stacking (see Fig. [2jc))PH These 
differences in the 2D band feature are best captured by 
the Lorentzian components of their fits. In particular, 
the Lorentzians with the highest intensities - i.e., cen- 
tered around 2685 cm -1 (green) and 2715 cm -1 (purple) 
- have very different intensities in the ABC samples (TGI 
and TG2), whereas they have almost equal intensities in 
the ABA trilayers (TG3), in agreement with the obser- 
vations reported in Ref. [2T] , 

We now turn our attention to the transport properties 
of trilayer graphene in double-gated transistor structures. 
This device geometry allows the independent control of 
the Fermi energy and the external perpendicular electric 
field E ex applied to the trilayers. In particular, the E ex 
is given by E ex = V tg /d tg — Vb g /d hg with V tg and Vbg 
the top- and back-gate voltages respectively, and d tg and 
dbg the thicknesses of the top- and back-gate dielectric. 
Fig. |3ja) and (c) show the 2-terminal square resistance 
(R sc jof trilayers with different stacking orders, measured 
for fixed values of Vtg as a function of Vb g . In all cases 
R sq displays a maximum (i?™ ax ) corresponding to the 
charge neutrality in the system. Clearly, the evolution of 
^sq ax w ^h E ex is markedly different for the two stacking 
orders. For ABC trilayer i?™ ax increases with increasing 



£? ex , whereas the opposite behavior is observed for ABA 
trilayer, i.e. i?™ ax decreases with increasing E ex . In both 
cases the position in Vbg of shifts linearly with Vt g , 
reflecting the changes in charge density induced by the 
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FIG. 3: (Color online) Square resistance R sq as a function of 
back-gate voltage for different fixed values of top-gate voltage 
at 4.2 K, shown for (a) ABC trilayer and (c) ABA trilayer with 
thicknesses of top-gate dielectric of 90 and 15 nm, respectively. 
Figures on the right present the schematic band structure of 
(b) ABC and (d) ABA trilayer graphene (considering only 70 
and 71). Application of a perpendicular electric field opens 
an energy gap for the ABC trilayer, while it results in a band 
overlap for the ABA trilayer. Solid and dashed lines are with 
and without the external electric field, respectively. 



two gates. 

These results can be understood by the effect of the 
perpendicular electric fields on the band structure of 
ABA and ABC graphene trilayers. Theory predicts that 
the interlayer asymmetry induced by the electric field 
opens an energy gap in the band structure of ABC tri- 
layers (Fig. gb))pMiIl whereas it causes a band over- 
lap in ABA trilayers (Fig. |3|d)). 13 The electric-field tun- 
able energy dispersion is a unique characteristic of few- 
layer graphene materials, and it paves the way to devices 
with unprecedented functionalities. Recent experiments 
in double-gated bilayer transistors have demonstrated an 
on/off current ratio of 100 at room temperature.^ On the 
other hand, very little is known experimentally on the 
electric-field tunability of the band structure of thicker 
few-layers and their stacking dependence. 

Fig. [4] illustrates the effect of a perpendicular magnetic 
field on the transport properties of ABA and ABC tri- 
layer graphene at T = 4.2 K. For ABC-stacked trilayer 
(TGI, fi « 1900 cm 2 V _1 s _1 ) the 2-terminal magnetore- 
sistance shows well-defined QHE plateaus at v = 6 and 10 
for B > 20 T (Fig.ga)). The filling factor v = n s </> £ _1 , 
where 0o is the flux quantum, matches well with the car- 
rier density n s = a (Vbg — ^cnp) calculated using a = 
7.2 x 10 10 cm _2 V _1 . The observed plateaus are expected 
from the 3-fold degenerate zero-energy LLs of the ABC 



trilayer graphene (E n ex B 3 / 2 ^n(n — l)(n — 2)) with 4- 
fold spin and valley degeneracy. We find QHE plateaus 
only away from the charge neutrality point (CNP) lo- 
cated at Vcnp ~ 20 V for V tg = (Fig. gb)). An ad- 
ditional plateau develops at v = 12, rather than at the 
expected v = 14, upon further increase of E ex (for exam- 
ple at Vb g = -90 V with V tg = -5 V). This observation 
suggests lifting of the valley degeneracy induced by the 
interlayer potential asymmetry,^ imposed by the top and 
back gates. 

The ABA-stacked trilayer device (TG3, 
fi w 1100 cm 2 V _1 s _1 ) develops QHE plateaus at 
v = 2, 4, 6, and 8 with a step of Ais = 2 (Fig. gc)). 
This observation is consistent with a recent theoretical 
prediction which includes the complete set of hopping 
parameters up to the next-nearest layer couplings 72 and 
75.^ This extended model predicts relative energy shifts 
of the monolayer- like and bilayer-like LLs in the ABA 
trilayer and a valley split of the zero-energy LLs by the 
72 and 75. As a result, the 12-fold zero-energy levels (4 
and 8 zero-energy levels from the monolayer-like and the 
bilayer-like subbands, respectively) split into 6 different 
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FIG. 4: (Color online) (a) Two-terminal magnetoresistance 
of the ABC-stacked trilayer at 4.2 K, shown for various back- 
gate voltages at Vt g = 0, except for a trace indicated. QHE 
plateaus develop at v — 6, 10, and 12 (dashed lines). The 
small deviations from the dashed lines reflect the contact re- 
sistance of our device (TGI), (b) Conductance G of the ABC 
trilayer device (TGI) as a function of Vb g (or carrier density 
n s ) at Vtg = 0. The mobility, /x, is estimated from the linear 
Vbg dependence of G at large back-gate voltages, (c) Magne- 
toresistance of ABA-stacked trilayer (TG3) shown for three 
different carrier densities at 4.2 K. (d) Schematic low-energy 
band structure (left) and zero-energy LLs (right) of ABA tri- 
layer graphene with all hopping parameters 70 — 75 included. 
The next-nearest layer couplings 72 and 75 shift the energy 
of monolayer-like (M) and bilayer-like (B) bands relative to 
each other, and also split zero-energy LLs into valleys. 
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energies with twofold spin degeneracy (Fig. |2Jd)), 
leading to the QHE plateaus at filling factor intervals 
of Az^ = 2. In addition, the presence of the external 
electric field generally splits the valley degeneracy of 
the LLs by the induced interlayer asymmetry.^ As 
opposed to the case of ABA trilayer, the electric-field- 
induced valley splitting is expected to be smaller for 
the inversion-symmetric ABC trilayer. Therefore, the 
4-fold spin and valley degeneracy is retained for the 
ABC trilayer device, resulting in QHE plateaus at v = 6 
and 10. Under the large external electric field, however, 
the valley splitting leads to the QHE plateau at v — 12. 

During the writing of our manuscript, we became 
aware of preprints dealing with QHE in ABA^ and 
ABC trilayersPHU] In Ref. [24], QHE plateaus in the 
ABA trilayer are observed at v = ±2, ±4, -6 but not at 
v = +6. The absence of a plateau at v = +6 is at- 
tributed to LL crossing. As pointed out in the Ref. [24] . 
actual plateaus developing in the ABA trilayer can de- 
pend on a specific B (or in our case Vbg) where mea- 
surements are performed, due to the LL crossing between 
the monolayer-like and the bilayer-like subbands. For the 
ABC trilayer, Refs. [25] and [27] report QHE plateaus at 
v = ±6, ±10, ±14..., consistent with our results except 



the plateau at v — 12, whereas Ref. [26] observed rather 
unexpected plateaus at v = ±9, ±18, and -30. 

In summary, we have investigated transport properties 
of trilayer graphene with different stacking order. Sam- 
ples with ABA and ABC stacking differ characteristically 
in the sequence of quantum Hall plateaus, in agreement 
with recent theory. The stacking order provides an addi- 
tional degree of freedom to tune the electronic properties 
of trilayer graphene, combined with the interlayer asym- 
metry controlled by top and back gates. 
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